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Abstract. Based on a power spectrum analysis of the IRAS ISSA maps, we present the first detection of the 
Cosmic far-Infrared Background (CIB) fluctuations at 60 and 100 /im. The power spectrum of 12 low cirrus 
emission regions is characterized by a power excess at spatial frequencies higher than k ~ 0.02 arcmin -1 . Most of 
this excess is due to noise and to nearby point sources with a flux stronger than 1 Jy. But we show that when these 
contributions are carefully removed, there is still a power excess that is the signature of the CIB fluctuations. The 
power spectrum of the CIB at 60 and 100 fim is compatible with a Poissonian distribution, at spatial frequencies 
between 0.025 and 0.2 arcmin -1 . The fluctuation level is ~ 1.6 x 10 3 Jy 2 /sr and ~ 5.8 x 10 3 Jy 2 /sr at 60 and 100 
/xm respectively. The levels of the fluctuations are used in a larger framework, with other observationnal data, 
to constrain the evolution of IR galaxies (Lagache et al. 2002). The detections reported here, coupled with the 
level of the fluctuations at 170 /im, give strong constraints on the evolution of the IR luminosity function. The 
combined results at 60, 100 and 170 /im for the CIB and its fluctuations allows, on the CIB at 60 /im, to put a 
firm upper limit of 0.27 MJy/sr and to give an estimate of 0.18 MJy/sr. 
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1. Introduction 

The Cosmic Far-Infrared Background (CIB) was detected 
by Puget et al. (1996), exploiting COBE-FIRAS data, and 
has now been firmly established over a range of wave- 
lengths (e.g. Dwek et al. (1998); Gispert et al. (2000); 
Hauser & Dwek (2001)). The intensity is quite high with 
respect to predictions based on evolutionary models of 
star formation in galaxy populations inferred from optical 
data. Source counts obtained with SCUBA at 850 /im and 
ISO at 170 and 15 /jm have partly resolved the CIB and 
shown a strong cosmological evolution. In the future, far- 
IR and sub-millimeter telescopes from ground and space 
will perform deep surveys over small areas, aimed at re- 
solving a substantial fraction of the CIB and to shed light 
on the number density, luminosity and spectral evolution 
of the infrared galaxy populations. However, investiga- 
tion of the clustering of these populations requires sur- 
veys over much larger areas. One way to tackle the limi- 
tation on the number of detected galaxies per field is to 
search for CIB fluctuations. So far, CIB fluctuations has 
only been observed at 170 /im in the FIRBACK fields 
(Lagache & Puget 2000) and at 90 and 170 /mi in the 
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Lockman hole (Matsuhara et al. 2000) surveys. These de- 
tections were probably dominated by the Poissonian con- 
tribution. By analyzing larger FIRBACK fields, power 
spectra seem to reveal correlated fluctuations, well above 
the cirrus contribution (Puget & Lagache 2002). These 
results are currently under further investigation. Above 
170 /im, only SCUBA observations might currently be able 
to give information on the fluctuations. However, no signif- 
icant CIB correlations have been detected in the SCUBA 
maps (Peacock et al. 2000). The CIB anisotropics are 
mainly contributed by moderate to high redshift star- 
forming galaxies, whose clustering properties and evolu- 
tionary histories are currently unknown. Since the clus- 
tering strength depends on the bias at the relevant red- 
shift, observing the CIB correlated fluctuations will pro- 
vide valuable informations on bulge and elliptical forma- 
tion, as well as potentially QSOs, thereby providing clues 
on the physical relations between dark matter and star- 
burst galaxies. 

In this paper we present the results of a power spectrum 
analysis of the IRAS 60 and 100 /im emission of 12 regions 
in the sky with very low interstellar emission. All these re- 
gions are characterized by an excess of power at high spa- 
tial frequencies with respect to the interstellar emission. 
We will show that this excess is not of instrumental origin 
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and may be attributed to the CIB. The discovery of the 
CIB fluctuations at 60 and 100 /im could give constraints 
on the number counts below the IRAS point source detec- 
tion limit. However, instead of extrapolating the number 
counts at 60 and 100 /im using the level of the fluctua- 
tions, which is of limited cosmological interest, we prefer 
to use the two detections in a more general framework of 
the modelisation of the IR galaxy evolution that combines 
all existing number counts, redshift distributions, and ob- 
servations of the CIB and its fluctuations, in the whole IR 
and submm range Lagache et al. (2002). 
The paper is organized as follow. After a presentation^ |J) 
of the data used for this analysis we describe in § || the 
power spectrum at 60 and 100 /im of the 12 fields selected. 
In § H we present the method used to separate the Galactic 
and extra-Galactic contributions to the power spectrum 
and we discuss our results in § |^. 

2. Data - the ISSA maps 

2.1. The IRAS Sky Survey Atlas 

The IRAS mission was originally designed to measure 
point sources but it also provided spectacular images of 
the Galactic dust diffuse emission. These images, consti- 
tuted of 500 x 500 pixels with a pixel size of 1.5' (i.e. the 
size of the field is 12.5° x 12.5°), are gathered in the IRAS 
Sky Survey Atlas (ISSA). Each ISSA map is the result 
of the combination of up to three individual maps (named 
for HCON Hours CONfirmation) . The HCON images were 
constructed from different observations of the same region 
separated by several months. The angular resolution of the 
60 and 100 /im band images are respectively 1.5' x 4.7' and 
3.0' x 5.0'. To study the power spectrum of the CIB emis- 
sion at 60 and 100 /im, we have worked on 12 ISSA maps 
with particularly low cirrus emission. 

2.2. Calibration of the ISSA maps 

Special care has been taken to have a consistent diffuse 
emission calibration through the Atlas. But, as the IRAS 
mission was designed to provide absolute photometry only 
for point sources, the ISSA images give only relative pho- 
tometry and cannot be used to determine the absolute 
surface brightness for diffuse emission. Based on a compar- 
ison with the DIRBE data, it has been shown (Wheelock 
et al. 1993; Schlegel et al. 1998) that, at large scale (> a 
few degrees), the amplitude of the fluctuations are overes- 
timated in the ISSA maps by a factor 1.15 at 60 /im and 
1.39 at 100 /im. We have thus applied these factors to our 
maps prior to the analysis. 

There is also an uncertainty on the zero level of the 
ISSA maps, for which the zodiacal emission has been sub- 
tracted (Wheelock et al. 1993). This uncertainty is dom- 
inated at 60 and 100 /im by an imperfect knowledge of 
the detector offsets and of the zodiacal emission. As we 
are looking at fluctuations of the signal, a global additive 
offset has no impact on the power spectrum. On the other 
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Table 1. Central pixel coordinates (ecliptic and Galactic) 
of the 12 ISSA maps selected for our analysis. 



hand, an imperfect zodiacal emission correction applied 
to the ISSA maps will have an impact on the large scale 
structure of the maps and then on the power spectrum 
(for low k values) . To restore appropriately the large scale 
structure of the selected fields we have compared each 
ISSA map with the DIRBE data (for which the zodia- 
cal emission correction was better done). The ISSA maps, 
multiplied by the appropriate gain value (0.87 at 60 /im 
and 0.72 at 100 /im) were convolved by the DIRBE beam 
and then subtracted from the DIRBE data. This offset 
map is then added to the ISSA map. Note that we do not 
use the Schlegel et al. (1998) IRAS rescaled maps, con- 
structed in a similar manner, for which it is impossible to 
recover the associated instrumental noise. 

2.3. Fields with low cirrus emission 

In this paper we present a power spectrum analysis of 
twelve high latitude fields. These fields were selected on 
the basis of their low cirrus emission (mean brightness 
~ 1 MJy/sr at 100 /im) but also on their redundancy; we 
have selected only fields for which each sky position has 
been observed at least twice in order to be able to esti- 
mate the contribution of the noise to the power spectrum. 
The typical stripping of the ISSA maps is very efficiently 
removed from the signal power spectrum by the noise esti- 
mate procedure using the difference between observations 
of the different HCON as will be shown in § 3.3 and Fig. 
2. 

The 60 and 100 /im ISSA maps of the twelve selected fields, 
gain and offset corrected, are shown in Figs |A.2| to [A.3| 
and their central coordinates are gathered in Table [l]. Five 
fields are located in the southern hemisphere, most of them 
in the neighborhood of the Marano field (Marano et al. 
1988), and seven are spread over the northern hemisphere 
(ISSA map number 376 contains the Lockman Hole). 



3. Power spectrum of the 60 and 100 /im sky 
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3.1. Power spectrum computation 



Therefore, the noise level of the ISSA map is given by: 



For a given image, the power spectrum P(k) is the ab- 
solute value of the Fourier Transform (we use the Fast 
Fourier Transform (FFT) function of IDL) of the image, 

averaged over constant values of k = 



k-x + ky. To make 
sure that map boundaries are not contaminating the com- 
puted power spectrum, we apodize the image with a co- 
sine function before computing the Fourier Transform. 
This prevents the introduction of strong discontinuities 
when the image is periodized in the FFT algorithm (see 
Appendix ^) . 

The power spectrum of the 100 (im ISSA map of one 
of our fields is shown in Fig. ^. Since the zodiacal emission 
has been subtracted in the ISSA maps, there are four main 
contributions to the power spectrum at these wavelengths: 
the cirrus emission, point sources, the CIB (resolved and 
unresolved) and the noise. Therefore, and if the noise and 
the signal are not correlated, the power spectrum P(k) 
can be expressed in the following manner: 

P(k) = 7(fc) [Pdurt(fc) + Source (k) + Pcib(fc)] + N(k) (1) 

where Pdust(fc), P c ib(k), P S ourcc(k) and N(k) are respec- 
tively the power spectrum of the dust emission, of the CIB, 
of individually detected point sources and of the noise. The 
factor j(k) represents the instrumental function. 



^diff 

vW2' 



(4) 



For the example given in Fig. [l], the noise level is (Ji ssa = 
0.048 MJy/sr. 

In general three HCON maps were averaged but parts 
of them may be undefined. Therefore, for each ISSA map 
we looked at the three corresponding HCON maps and 
built a mask n(x,y) that gives the number of defined val- 
ues that were averaged at each positions. Then, the noise 
map of a given ISSA map is estimated by subtracting two 



HCON maps and divide the result by y/2Wn(x,y). 

The power spectra of the ISSA map and of its noise 
map are shown in Fig. |l](bottom-left). Both power spectra 
meet at small scales (k ~ 0.2 pixel -1 ) where the signal is 
noise dominated. The noise map is characterized by a ~ 
£-0.8 power spectrum. We have investigated the variations 
of the statistical properties of the noise on ISSA maps 
where three complete HCONs exists. By looking at the 
three possible difference maps built from three HCONs, 
we found that the shape of the power spectrum of the 
noise is rather constant with time and that the absolute 
level of the noise varies by less than 15% from one HCON 
to the other. This is equally true at 60 and 100 /im. 



3.2. Noise power spectrum 

One of the main limitation of the component separatio 
to the power spectrum is the estimate of the noise cor 
tribution. For ISSA maps, this contribution can be est 
mated accurately as they are the combination of up t 
three HCON maps. Indeed the noise power spectrum c 
the HCON maps, and therefore of the ISSA map, can b 
estimated by subtracting two HCON maps of the sam 
region. An example of the difference between two HCO] 
maps for a typical low brightness ISSA map is shown i 
Fig. [I] (top-right). The characteristic stripping of the IRA 
data is seen in this difference map. The distribution of di: 
ference values is symmetric (when there is no strong poir 
sources) and well fitted by a sum of two Gaussian func 
tions (see Fig. [j], bottom-right). 

As the distribution of the difference values is the sui 
of two Gaussian functions with an equivalent width <7di 
and assuming that the noise in the IRAS survey is sta- 
tionary to a good approximation (as tested below), we 
can conclude that the noise of an individual HCON can 
also be characterized by an equivalent width: 

Cdiff /r,\ 

°hcon — — 7=-- ( z ) 

Furthermore, the noise of the ISSA map, built from n 
HCON maps, is also Gaussian distributed with a width 
of: 
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Fig. 2. Power spectrum (black solid line) of a typical 
bright ISSA field with no strong point sources. The cirrus 
contribution to the power spectrum has been fitted by 
a power law (P(k) cx k~ 2m for k < 0.015 arcmin -1 - 
dash line) . The dotted line is the cirrus fit multiplied by a 
Gaussian PSF of a — 0.07 arcmin -1 . The grey line is the 
power spectrum of the ISSA map (black solid line) divided 
by the Gaussian PSF. 
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(3) 
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Fig. 1. Top left: The 100 /im ISSA map number 47. Top right: Noise map of the 100 /im map on the left. This noise 
map has been obtained by subtracting two individual HCON images that were used to build the ISSA map (see text 
for details). Bottom left: Power spectrum of the ISSA map (upper curve) and of the noise map. The power spectrum 
of the noise map has been fitted by a power law oc fc~ - 83 for scales k < 0.1 arcmin" 1 . Above k=0.1 arcmin -1 , the 



drop of the power spectrum is due to the instrumental function (see § |3.3| ). Bottom right: Histogram of the noise 
values, fitted by a sum of two Gaussian functions (Amplitudei = 1376, o\ — 0.04 MJy/sr; Amplitude2 = 711, a<i = 
0.06 MJy/sr). 



3.3. Point spread function 

The multiplicative factor ^{k) in Eq. [j] represents the effec- 
tive instrumental filter function of the system (telescope, 
focal plane assembly, sky scanning and map making). For 
some ISSA maps the instrumental function can be esti- 
mated directly from the power spectrum of the map. This 
is the case for maps where the brightness fluctuations are 
dominated by the cirrus emission and where the power 
spectrum is simply: 



Following Gautier et al. (1992), the power spectrum of the 
cirrus dust emission follows a power law (Pdust(^) = Ak^). 
This behavior has been confirmed through a detailed 
study of Hi 21 cm emission by Miville-Deschenes et al. 
(2002). These studies show that although the power spec- 
trum slope may vary from region to region by ~20 % there 
is no systematic break at high spatial frequencies. On the 
basis of these studies, we thus conclude that the power 
law representation for the cirrus power spectrum does not 
introduce any systematic effect. 



P(k) = 7 (fc) x P dust (k). 



(5) 
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To estimate j(k) directly on the power spectrum of 
ISSA maps, we have selected 20 relatively bright regions 
(mean brightness greater than 8 MJy/sr) where the noise 
and CIB contributions are negligible. To make sure that 
the fluctuations are dominated by the cirrus emission, 
point sources were filtered out using a median filtering. 
An example of the power spectrum of such an ISSA map 
is shown in Fig. |^. We found that the effective IRAS in- 
strumental filter function in Fourier space is well describe 
by a Gaussian function both at 60 and 100 /zm. To de- 
termine the width of the instrumental function, the power 
spectrum of the 20 fields were fitted by the following equa- 
tion: 



P(k) = exp 



2*2 



Ah? 



(6) 



which represents the cirrus power law multiplied by the 
Gaussian instrumental function. We have found <Jk — 
0.073 ± 0.005 at 60 fan and a k = 0.065 ± 0.005 at 100 
fim. In the real space, this corresponds to a Gaussian in- 
strumental function of a — 1.5 ± 0.1 arcmin at 60 //m and 
a = 1.7 ± 0.1 arcmin at 100 /im. [j] 

The cutoff seen in Fig. || at frequencies larger than 
0.05 arcmin -1 is where the effect of the instrumental filter 
function is expected. When the power spectrum is divided 
by the Gaussian instrumental function (grey line) one sees 
that this modeling does not stand for scales smaller than 
4.5 arcminutes (k > 0.22 arcmin -1 ). This is due to the 
effective resolution of the IRAS data. The power spectrum 
has to be cut at the Nyquist frequency corresponding to 
the 1.5' pixel ISSA map. However, the effective resolution 
of the IRAS data is about 3.8' and 4.25' respectively at 
60 and 100 /mi. For that reason, all the power spectrum 
at frequency above 0.1 arcmin -1 , which corresponds to 
the Nyquist frequency of the effective resolution, has to 
be taken with cautious. 

4. Separation of the Galactic and extra-Galactic 
components 



In Figs. A. 2 to A. 2 we present the power spectrum of the 



60 and 100 /im ISSA maps of our twelve fields (black line). 
For all these fields, the power spectrum of the noise has 
been subtracted and the result has been divided by the 
PSF 7(fc). At large scale, all the power spectra follow a 
power law typical of cirrus emission (Gautier et al. 1992). 
But at k ~ 0.01 arcmin -1 we notice a break in the spectra 
and an excess from the cirrus power law at small scales. 
This power excess is observed for all our fields, at 60 and 
100 fim. 

4.1. Contribution of strong point sources 

The power excess observed at small scales reveals the 
presence of a component with a flatter power spectrum, 



1 This is compatible with what we found by fitting a 
Gaussian function on point sources. 



like what would be expected from noise or randomly dis- 
tributed point sources. On the one hand, at this stage, 
a noise contribution is unlikely as it has been well esti- 
mated and removed accordingly. On the other hand, in 
fields with such low cirrus emission, several extra-galactic 
point sources were detected. In each of our fields, between 
200 and 300 point sources are listed in the IRAS Point 
Source Catalog (PSC), the vast majority of them being 
extra-galactic objects. We remove the contribution of all 
point sources with flux greater than 1 Jy at 100 /im (and 
the same sources at 60 /mi). The cut at 1 Jy is chosen such 
that the PSC above this limit will be nearly complete. This 
gives a well define separation between the sources con- 
tributing to the unresolved CIB and the resolved sources. 
This is necessary for any quantitative use of our results. 
We choose this cut in flux as low as possible to remove the 
contribution of the nearby and resolved galaxies that dom- 
inate the power spectrum at scales k > 0.01 arcmin -1 (for 
those galaxies, the spatial distribution is indistinguishable 
from a Poissonian distribution). Moreover, according to 
recent simulations, the cut of 1 Jy at 100 /im almost corre- 
sponds to the flux where the non-euclidean part appears in 
the number counts. To remove the contribution of strong 
point sources (Iioo^m > 1 Jy) to the power spectrum, they 
were filtered out from the ISSA map. For each ISSA map, 
we have extracted the point sources of the IRAS Point 
Source Catalog with a flux greater or equal to 1 Jy. Then 
we have applied a median filtering to the ISSA map to 
estimate the background level and, at each point source 
position, we have removed the points in a 12 x 12 window 
that were more than 4 times the noise level (cissa) above 
the background. Finally, the missing points were replaced 
by a bilinear interpolation. 

The power s pect ru m of the filtered ISSA maps is the 
blue line in Figs A. 2 to A. 3. A power excess at small scales 
(k > 0.02 arcmin -1 ) is still apparent on all power spectra, 
at 60 and 100 /im. This excess is likely to be the signa- 
ture of the unresolved cosmic far-infrared background. To 
characterize it, the contribution of the cirrus emission has 
to be removed. 



4.2. Contribution of the cirrus emission 

As we do not have any independent tracers of the cirrus 
emission at the IRAS scales for the whole fields, the only 
way to derive the cirrus contribution is to use the sta- 
tistical properties of their spatial distribution. Following 
Gautier et al. (1992) and Miville-Deschenes et al. (2002) 
the power spectrum of interstellar dust emission is well de- 
scribe by a power law. To remove this contribution to the 
power spectrum of the filtered ISSA maps, we have fitted 
a power law on the large scale part of the power spec- 
trum of the filtered maps. The best compromise between 
the number of point on which to fit the power law and 
the contamination from the CIB was to fit on k < 0.02 
arcmin -1 . The point corresponding to the largest scale 
was not used for the fitting as it suffers from a large sta- 
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tistical error. The result of the fitting, and its associated 
uncertainty, are shown for each 1SSA map in Figs |A.2| to 
. The uncertainty on the slope reflects the dispersion 
of the points around the fit. The red line in these figures is 
the result of the subtraction of the cirrus power law from 
the power spectrum of the filtered map (blue line); it is 
the signature of the CIB. 

4.3. Contribution of the cosmic far-infrared 
background 

The average of the 12 CIB contributions is shown in Fig. || 
(black lines). The blue lines in this figure indicate the la- 
uncertainty on the determination, including the PSF and 
the estimate of the cirrus contribution errors. The la un- 
certainty is very close to the rms variations of the 12 av- 
eraged values. The dispersion between the fields is thus 
dominated by the PSF and cirrus contribution errors, we 
do not see any significant variation from field to field. 
Moreover, it is important to note that the detected excess 
at small scales is uncorrelated with the cirrus emission (the 
excess disappears when the cirrus contribution increases). 
Therefore, our fields being spread over the sky, the excess 
at small scale is compatible with isotropic properties. This 
is an other argument in favor of an extra-galactic origin 
for the excess. 

We can note on Fig. A. 2-9 that the power spectrum of 
the resolved IRAS sources is varying by factor up to 10 
from field to field. The CIB (resolved and unresolved) fluc- 
tuations are dominated by the brightest sources in ISSA 
maps. The number of such sources is quite small, leading 
to a high cosmic variance, and thus large variation from 
field to field. On the contrary, for the unresolved CIB fluc- 
tuations, as sources above 1 Jy (about 300 per field) are 
removed, the cosmic variance is lower than 6%, which is 
negligible with respect to the la uncertainty. 
The CIB is detected over one decade in scales, from 5 to 
50 arcminutes. The power spectrum of the CIB at 60 and 
100 /xm is compatible with a Poissonian distribution at 
spatial frequencies between 0.025 and 0.2 arcmin— 1. The 
fluctuation level is ~ 1.6 x 10 3 Jy 2 /sr and ~ 5.8 x 10 3 
Jy 2 /sr at 60 and 100 respectively. 

5. Discussion 

5.1. CIB fluctuations at 60 and 100 fim 

We have shown here that the power spectrum of high lati- 
tude fields at 60 and 100 (im is characterized by a break at 
small scales (near 0.02 arcmin -1 ). As indicated by Gautier 
et al. (1992), the power spectrum of the dust emission 
is usually well described by a power law proportional to 
£,-3±o.2_ ^ m ore detailed study of the statistical proper- 
ties of the interstellar cirrus HI 21 cm emission has been 
carried out by Miville-Deschenes (1999). In this analysis 
it is shown that there are limited variations of the spec- 
tral index of the power law from field to field but, what 
is most important for the present work, the power spec- 



trum of cirrus emission for scales smaller than 12.5° is 
always characterized by a single power law with no break. 
Therefore, it is unlikely that the power excess observed 
here at small scales is of interstellar origin. 

We could also wonder if this break is of instrumental 
origin. It was shown by Wheelock et al. (1993) that the 
response of the IRAS detectors are affected by memory ef- 
fects. This produces variations of the detector response as 
a function of scale. This effect is more important at small 
scales (under a few tens of arcminutes) but Wheelock et al. 
(1993) have shown that the amplitude of the fluctuations 
at these scales were underestimated. This effect will thus 
produce a drop of the power spectrum at small scale and 
cannot explain the power excess observed here. 

On the other hand, instrumental noise could produce 
such an excess in the power spectrum. But, as the IRAS 
ISSA maps result from redundant individual observations, 
we were able to estimate the contribution of the noise to 
the power spectrum. We are aware that our estimate of 
the noise rely on the fact that the individual HCONs are 
independant. This is not perfectly true as, in the construc- 
tion of the HCONs, an offset was added to each scan to 
minimize the difference between different observations of 
the same position. Therefore the noise level estimated by 
subtracting HCONs may be underestimated at the scale 
of a scan, which is a few degrees. But at this angular scale 
the signal is completely dominated by the cirrus emission, 
even in the low brightness regions selected for our anal- 
ysis. At the scale of a few arcminutes where the CIB is 
detected, the noise contribution to the power spectrum 
has been removed accurately. 

In fact it appears that most of the power excess can be 
attributed to the numerous extra-galactic point sources 
that are present in such a low cirrus emission field. When 
the strong {Iwo^m > 1 Jy) point sources are removed 
from the ISSA maps, we recover a power spectrum 
typical of cirrus emission at low spatial frequencies 
but with still a power excess at small scales (k> 0.02 
arcmin" 1 ) that can be attributed to the unresolved 
cosmic infrared background. Moreover, the residue has 
homogeneous properties over the sky, consistent with CIB. 

Knox et al. (2001) computed the expected power spec- 
trum of the CIB at several frequencies (v < 1060 GHz), 
exploiting the far-IR volume emissivity derived from the 
count models of Guiderdoni et al. (1998) and assuming a 
bias 6 = 3, constant with redshift. They concluded that 
the clustering-induced fluctuations can match those of the 
CMB at I <300. They also predict a shape of the CIB 
power spectrum peaking around scales of 1-3 degrees. This 
broad maximum, if present, is at the limit of our frequency 
range where the noise is becoming large, making the de- 
tection of the clustering very difficult. The power spec- 
tra of the CIB at 60 and 100 /zm are compatible with a 
Poissonian distribution with levels ~ 1.6 x 10 3 Jy 2 /sr and 
~ 5.8 x 10 3 Jy 2 /sr respectively. 
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Fig. 3. CIB contributions to the 60 (left) and 100 /zm (right) power spectra. The black line is the average of the CIB 
contributions for the 12 fields. The blue lines indicates the uncertainty on the CIB contribution, computed from the 
uncertainty on the cirrus slope and on the PSF width. 
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Fig. 4. Redshift distribution of the sources making the CIB and the fluctuations at 60, 100 and 170 /im. The first 
two-left panels are for the total contributions, the last two-right panels, for sources with flux below 1 Jy (Lagachc 
et al. 2002) 
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5.2. CIB intensity and anisotropy amplitudes color 
ratio 

The CIB rms fluctuations in the IRAS maps corresponding 
to the white noise power spectra are: 

a 2 = J P cib (k)2nkdk Jy 2 /sr 2 (7) 

giving cr=0.048 MJy/sr and ct=0.09 MJy/sr at 60 and 
100 /jm respectively. As shown by Gispert et al. (2000), 
the CIB at different wavelengths is dominated by sources 
at different redshifts, larger wavelengths being dominated 
by more distant sources. The same applies to the fluctua- 
tions: lower frequencies probe higher redshifts (e.g. Knox 
et al. (2001)) The fluctuations at 60 /im are dominated by 
nearby bright objects. When we remove these objects, the 
residual fluctuations are quite low. At 100 /im, the contri- 
bution of higher redshift objects is increasing, leading to 
higher level of residual fluctuations. Therefore, the ratio 
of the 60 to 100 fluctuation is qualitatively consistent 
with what is expected. 

This is illustrated more quantitatively in Fig.^ on 
the panels showing the redshift distribution of sources 
contributing to the CIB intensity and the fluctuations 
(Lagache et al. 2002). It is clear from these figures that 
the z-distribution of the fluctuations is bimodal, with one 
contribution at redshift lower than 0.25 and the other one 
centered at redshift around 1. For all sources the ratio 
of nearby to moderate-redshift source contribution to the 
fluctuations is equal to 3.8, 1.8 and 1.3 at 60, 100, 170 
/zm respectively, illustrating that fluctuations at larger 
wavelengths are dominated by more distant sources. When 
the brightest sources are removed the ratio of nearby to 
moderate-redshift contribution becomes equal to 1.7, 1 
and 0.74 at 60, 100, 170 /im. In this case, at 100 /jm the 
contribution to the fluctuations of nearby and moderate- 
redshift sources is the same, becoming lower at higher 
wavelength. At 60 /im, the fluctuations are still dominated 
by the nearby objects. 

For the three wavelengths, the CIB is mainly due to 
sources at redshift around 1. A detailed analysis of the CIB 
fluctuations at 100 and 170 /im (which is beyond the scope 
of this paper) will give information on the distribution of 
sources at z^l making the bulk of the CIB. This is partic- 
ularly true at 170 /zm where sources with flux lower than 
4cr=135 mJy can be removed (Dole et al. 2001), leading 
fluctuations highly dominated by the moderate-z sources. 

We can compute the ratio of CIB fluctuations to in- 
tensity (R\ = j^) at 100 /im and compare it with the 
previous determination at 170 /zm. To compute Rioo and 
R170, we use: 

— a CIB intensity at 100 /im of 0.5 MJy/sr (Renault et al. 
2001), leading to Ri 00 = 0.18 

— a CIB intensity at 170 /im of 1 MJy/sr Lagache & Dole 
(2001). For the fluctuations, we assume that a cut of 
1 Jy at 60 and 100 /im corresponds to the same cut 
at 170 /im. We obtain rms fluctuations of 0.12 MJy/sr 



(corresponding to ^11000 Jy 2 /sr, (Puget & Lagache 
2002)). This gives Ri 70 = 0.12 

The ratio R is decreasing between 100 and 170 /im, as ex- 
pected. At 60 /im, the minimal hypothesis is to consider 
that Reo is equal to Rioo=0.18 which gives an upper limit 
of 0.27 MJy/sr on the CIB intensity. An illustrative log- 
arithmic extrapolation of R out to 60 /im gives R6o=0.27 
leading to an estimate of the CIB intensity at 60 /zm of 
0.18 MJy/sr, which is significantly smaller than the pre- 
vious determination of Finkbeiner et al. (2000) of 0.56 
MJy/sr. 

5.3. Implication for the component separation 

This work suggests that the high latitude IRAS maps, in 
the lowest cirrus regions, cannot be used as a tracer of the 
interstellar extinction structure as proposed by (Schlegel 
et al. 1998). In fact it should be noted that it is only above 
an intensity of order of 10 MJy/sr at 100 /im that the CIB 
fluctuations are lower than the cirrus contribution at the 
smallest scales (k~0.2 arcmin -1 ). 

Present and future CMB observations, above 100 
GHz, with high sensitivity bolometers need to remove 
foreground contributions (cirrus and CIB fluctuations). 
The CIB spectrum being significantly "colder" than 
the cirrus spectrum (/ C j rrus 100/Z77i// C i rrus lmm ~ 30; 
^ciBlOO/im/JciBlniw ~ 5), the relative contribution of 
the CIB will increase with wavelength. It is thus expected 
that at 1 mm the range in 1 space dominated by the CIB 
will be much more extended than sees at 100 /zm. This 
question will be dealt in a forthcoming paper. 
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Appendix A: Power spectrum computation 
method 

The power spectrum of an image f(x,y) of Fourier 
Transform f{k x ,k y ) is computed from the amplitude 
A(k x ,k y ) defined as 



A(k x , ky) — f(k x , ky)f {k Xl k y ) — \ f(k x , k y )\ 



(A.I) 



The power spectrum P(k) is the average of A(k x ,k y ) on 



annulus of constant k = * /fc.? 



k 2 

K y . 



The amplitude image 



of a typical ISSA map is shown in Fig A.l (top-right). 
In this figure, the frequency k — is a the center of the 
amplitude image. The amplitude decreases with increas- 
ing value of k but one also notices the presence of a cross 
that reveals an increase of the fluctuation level in the hor- 
izontal and vertical directions. This cross is caused by the 
Fast Fourier Transform algorithm that makes an "infinite 
pavement" with the image prior to compute the Fourier 
Transform. As the ISSA images are not periodic objects, 
this operation produces discontinuities where the left (top) 
meets the right (bottom) side of the image. These horizon- 
tal and vertical discontinuities are responsible for the cross 
seen in the amplitude image. The main problem with this 
cross is that it increases artificially the power when the 
average over k is done. Furthermore, the power increase is 
not constant as a function of k which modify the slope of 
the spectrum as well as the normalization. 

The usual method to get rid of this effect is to apodizc 
the image so that the left (top) and right (bottom) sides of 
the image have a similar flux level. Another method which 
do not modify the image but acts on the amplitude im- 
age itself is to compute P(k) by taking the median value 
at constant k. For ^85% of the ISSA maps we have in- 
spected, these two methods give very similar results. But 
when the cross is very strong, the power at small scales is 
still overestimated by the "median" method. Furthermore, 
at small k values, there are less points to average and a 
large fraction of them are affected by the cross effect; at 
these frequencies, the median value may not be significant. 
For all these reasons, we have adopted the "apodization" 
method. 
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Fig. A.l. Upper left: Typical ISSA map with a large 
scale gradient. Upper right: Amplitude of the Fourier 
Transform of the image on the left. The cross seen here 
is produced by the left-right and top-bottom flux differ- 
ence in the image. Lower right: Amplitude of the Fourier 
Transform of the image apodized with a cosine function 
to minimize the left-right and top-bottom discontinuities. 
Lower left: Power spectrum of the original image (solid 
line) and of the apodized image (red line) . The power spec- 
trum computed with the Median method is also shown 
(dotted line). It follows almost perfectly the power spec- 
trum of the apodized image. 
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Fig. A. 2. Power spectrum of a typical 60 (im ISSA maps. The black curve is the power spectrum of the raw ISSA 
map from which the power spectrum of the noise has been subtracted. The blue curve is the power spectrum (also 
noise subtracted) of the map from which point sources stronger than 1 Jy in the IRAS point sources catalog were 
removed. A power law fit of the blue curve (for k < 0.02 arcmin -1 ) is also shown. The red curve is the difference 
between the blue curve and the power law fit. 
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Fig. A. 3. Power spectrum of a typical 100 fxm ISSA maps. The black curve is the power spectrum of the raw ISSA 
map from which the power spectrum of the noise has been subtracted. The blue curve is the power spectrum (also 
noise subtracted) of the map from which point sources stronger than 1 Jy in the IRAS point sources catalog were 
removed. A power law fit of the blue curve (for k < 0.02 arcmin -1 ) is also shown. The red curve is the difference 
between the blue curve and the power law fit. 



